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We study properties of two-color QCD at imaginary chemical potential (/i) from the view point of 
the Roberge- Weiss (RW) periodicity, the charge conjugation and the pseudo-reality. At fi — inT/2, 
where T is temperature, the system has both the Z2 and the charge-conjugation symmetry. The two 
symmetries are preserved at lower T, but simultaneously broken at higher T. The Polyakov-loop 
extended Nambu-Jona-Lasinio (PNJL) model has the same properties as two-color QCD for the 
RW periodicity, the charge conjugation and the pseudo-reality. The nontrivial correlation between 
the chiral and Z2 symmetry breakings are investigated by introducing the entanglement vertex in 
the PNJL model. The order of the Z2 symmetry breaking at the RW endpoint is the second order 
when the correlation is weak, but becomes the first order when the correlation is strong. We also 
investigate the impact of the correlation on the phase diagram at real fi. 

PACS numbers: 11.30.Rd, 12.40.-y, 21.65. Qr, 25.75.Nq 



J2' 
9-1 

43 



> 

CO 
OO 
CN 
CN 

OO 
O 
CN 



X 



I. INTRODUCTION 

Elucidation of QCD at finite temperature (T) and fi- 
nite quark-number chemical potential (/i) is one of the 
most important subjects in hadron physics. Lattice QCD 
(LQCD) is the first-principle calculation, but has the 
sign problem at real fi. Particularly at fi/T > 1, the 
LQCD calculation is not feasible, although several meth- 
ods have been proposed so far to circumvent the problem; 
see for example Ref. [1]. For this reason, effective models 
such as the Polyakov-loop extended Nambu-Jona-Lasinio 
(PNJL) model [2-6] are widely used to investigate QCD 
at finite fi. 

Symmetries are important to understand QCD. In this 
paper we focus our discussion on global symmetries. 
QCD has the chiral symmetry in the limit of zero cur- 
rent quark mass (to) and the Zjv o symmetry in the limit 
of infinite m, where N c is the number of colors. The 
charge conjugation (C) symmetry is preserved at fi = 0, 
but not at finite fi. The QCD partition function Z(fi) is, 
however, invariant under the combination of the charge 
conjugation and the parameter transformation fi — > —fx. 
This leads to a relation between Z(fi) and Z(—fi): 



Z{fj) = Z(-fi). 



(1) 



The Zjv c symmetry is broken explicitly for light quarks 
of our interest, but it plays an important role at imagi- 
nary chemical potential ft, = i6T, where 9 is the dimen- 
sionless imaginary chemical potential. The QCD parti- 
tion function Z{9) is not invariant under the Zjv e trans- 
formation 



q 



Uq, A v -> UA U U~ 



i{d v U)U~ 



(2) 



where q is the quark field, A v is the gauge field and 
U (x, t) are elements of SU(7V C ) with the boundary con- 
dition U(x, f3 = \/T) = cxp(2i-Kk/N c )U{x, 0) for integers 



k from 1 to N c — 1. However, Z(9) is invariant under 
the combination of the Zjv c transformation and the pa- 
rameter transformation 9^9 + 2nk/N c , i.e., under the 
extended 1in c _ transformation [4] 



Uq, A v -> UA^U- 1 - i(d v U)U~ 



This means that 



Z(6) = Z{6 + 2nk/N c ). 



(3) 



(4) 



The relation (4) is called the Roberge- Weiss (RW) pe- 
riodicity [7]. QCD at finite is governed by the RW 
periodicity [4, 7]. 

The partition function Z{9) is £>-cvcn because of (1). 
This ^-evenness and the RW periodicity yield a constrain 
on 8 as < 9 < tt/N c and hence on fi 2 as 



l? > ~{Tir/N c ) 2 



(5) 



At the lower bound fi 2 = ~(Tir/N c ) 2 , i.e. at 6 = ir/N Cl 
the charge conjugation comes back to a symmetry. In 
fact, Z(ir/N c ) is changed into Z(—ir/N c ) by the charge 
conjugation, but Z(—tt/N c ) is identical with Z(tt/N c ) be- 
cause of the RW periodicity (4). This is true for any A?" c . 

The Zjv c transformation, meanwhile, becomes a sym- 
metry at 9 = tt/N c , only when N c = 2. In the case 
of N c = 2, the Z2 transformation changes Z(ir/2) into 
Z(7r/2±7r), but Z(tt/2±7t) is identical with Z(tt/2) ow- 
ing to the charge conjugation (1). This is not the case 
for N c > 3. For example in case of N c = 3, the Z3 trans- 
formation changes Z(ir/3) into Z(tt/3 — 2nk/3). The 
partition function Z(ir/3 — 2irk/3) is transformed back 
to Z(jr/3) by the charge conjugation when k = 1, but 
not when k = 2. At the lower bound fi 2 = —(Tn/N c ) 2 , 
QCD thus possesses both the C and the Zjv c symmetry 
when N c — 2, but has only the C symmetry when N c > 3. 
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At the lower bound /i 2 = —(Ttt/N c ) 2 , a first-order 
phase transition takes place for any N c [7]. The tran- 
sition is induced by the spontaneous breaking of the 
C symmetry [8]. The transition is now called the RW 
phase transition. A current topic on the RW phase tran- 
sition is the order at the endpoint. Recent three-color 
LQCD simulations show that the order is the first-order 
for small and large m, but the second-order for interme- 
diate m [9-12]. The order may be the first-order for both 
two-flavor [9, 10] and three-flavor cases [11, 12], when the 
pion mass has the physical value. If the order is the 
first-order, the RW endpoint becomes a triple-point at 
which three first-order transition lines meet. The PNJL 
model reproduces these results [13-15]. 

Two-color QCD has some interesting points. The num- 
ber of colors, N c , can vary from 2 to infinity, that is, 
realistic three-color QCD is between two-color QCD and 
large N c QCD. In this sense, understanding of both two- 
color and large N c QCD is important. The algebraic ap- 
proach based on the pseudo-reality [16] plays an impor- 
tant role in two-color QCD, while large N c QCD is well 
understood by the geometric approach based on the 1 /N c 
expansion [17] or the AdS/CFT correspondence [18]. In 
virtue of the pseudo-reality, two-color LQCD has no sign 
problem not only at imaginary /i but also at real /x [19], 
and consequently LQCD data are available there [20-25]. 
Furthermore, two-color LQCD has higher symmetries at 
fx = iTtt/2, as mentioned above. 

In this paper, we study properties of two-color QCD 
at imaginary \x from the view point of by the RW pe- 
riodicity, the charge conjugation and the pseudo-reality. 
Particularly, we clarify the relation between the Z 2 and C 
symmetries at 6 — ir/2. The PNJL model has the same 
properties as two-color QCD for the RW periodicity, the 
charge conjugation and the pseudo-reality. The PNJL 
model is then used to investigate two-color QCD con- 
cretely. Particularly, the nontrivial correlation between 
the chiral and Z2 symmetry breakings are investigated. 
In the model approach, we mainly consider the two-flavor 
(Nf = 2) case. 

This paper is organized as follows. In Sec. II, some 
properties of two-color QCD arc derived by using the 
RW periodicity, the charge conjugation and the pseudo- 
reality. In Sec. Ill, the two-color PNJL model is for- 
mulated with the mean-field approximation. Numerical 
results are shown in Sec. IV. Section V is devoted to 
summary. 



II. PROPERTIES OF TWO-COLOR QCD 

We consider the case of Nf = 1 for simplicity and 
assume that /i is either real or pure imaginary. The par- 
tition function Z of two-color QCD is obtained in Eu- 
clidean spacctime as 



with 



Z = J DAdf&[M{ji)]exp[~F%,] 



(6) 



M(fx) = D + m - 74/i, 



(7) 



where the Dirac operator D is defined by D = 7^(9^ — 
iA^) for the current quark mass m and the gauge field 
A v , and F^ u = d ll A u - d v A^ - i[A^, A v ], For later con- 
venience, we define Pauli matrices ti in color space and 
the Dirac charge-conjugation matrix C = 7274- 
The fermion determinant det[A/(/i)] satisfies 



(det[M(n)])* = dct[M (-/J,*)] 



(8) 



(det[M(//)])* = det[M(»t] = dct[7 5 M(/^ 75 ] 

= det[M (-//)]. (9) 

The relation (8) indicates that det[M(/x)] is real when /1 
is pure imaginary. The relation (8) is true for any N c . 
Two-color QCD has the pseudo- reality [16], 

Dt 2 C l5 =t 2 C l5 D\ (10) 

and hence the fermion determinant satisfies 

det[M(/i)] = dct[(t 2 C7 5 r 1 M(/i)(^7 5 )] 

= (det[M(//)])* (11) 

in virtue of the pseudo-reality (10). This means that 
det[M(/x)] is real when fj, is real. Two-color LQCD thus 
has no sign problem at both real and pure imaginary 
fi [19]. The charge-conjugation relation (1) is obtained 
from (8) and (11). 

The Polyakov loop $ is the vacuum expectation value 
of the Polyakov-loop operator 



tr c ( V exp 



l/T 



cLtAa 



(12) 



with the time-ordering operator V . The operator L is 
real, because 



V 



1 



tr c ( V exp 



l/T 



dr{A 4 y 



(13) 



where the second equality is obtained from the identity 
t 2 A„t 2 = -(A„) T = -(A,)*. 

The Polyakov loop $ is then real at both real and pure 
imaginary /i; note that det[M(//)] is real there. Under 
the charge conjugation, the factor — (A4) 1 - is transformed 
into A4, and hence the second equality of (13) means that 
L is C-invariant. Using this property and the relation (1), 
one can see that 



3(A0 = $(-p). 



(14) 



For pure imaginary chemical potential fi = iT8, it is 
convenient to introduce the modified Polyakov loop 



(0) = $(0)e 



w 



(15) 
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The modified Polyakov loop satisfies the RW periodicity 



*(0) = *(6» + tt) 



(16) 



since it is invariant under the extended Z^ transforma- 
tion [4]. It also satisfies 



#(0)* = 



(17) 



because of (14) and the realness of $. At = tt/N c = 
7r/2, both the Z 2 and C transformations become symme- 
tries, as mentioned in Sec. I. The order parameter of the 
C symmetry is a C-odd quantity. The relation (17) indi- 
cates that the imaginary part Im[4'(0)] is C-odd (0-odd) 
and hence an order parameter of the C symmetry [8]. 
Meanwhile, the order parameter of the Z 2 symmetry is 
$. At = 7r/2 the two order parameters are identical, 
since lm[^(0)] = $(0)sin(0) = $(0). This indicates that 
the spontaneous breakings of Z 2 and C symmetries take 
place simultaneously. 

Two-color LQCD simulations were made in Ref. [20, 
21, 23, 25] for the N f = 8 case. The LQCD results at 
= 7r/2 show that <t> = at small T but finite at large 
T. This indicates that the spontaneous and simultaneous 
breaking of Z2 and C symmetries occurs at some temper- 
ature T c . Further analyses are made in Sees. Ill and IV 
by using the PNJL model. 

For the Nf = 2 case, the fermion determinant is de- 
scribed by 



det[M(/x)] = det[M u {ii)]det[M d (n)}, 



(18) 



where det[M u (/x)] and dct [M^/i)] are the fermion deter- 
minants for u- and d-quark. Using the operator Ct 2 only 
for det[Md(/x)], one can get the relation 

det[M(n)) = det[M u ( M )]dct[(Ct 2 )^ 1 Md(M)(^2)] 

= det[M„(/i)]det[M d (-/i)]. (19) 

This relation indicates that the 1 + 1 system with finite fi 
is identical with the 1 + 1 * system with the same amount 
of isospin chemical potential Hi so . Hence, [i 2 so also has a 
lower bound 



(20) 



The diquark condensate in the former system corre- 
sponds to the pion condensate in the latter system. Be- 
cause of this symmetry, we consider the former system 
only in the present paper. 



III. PNJL MODEL 

We consider the N c = N f = 2 case. The PNJL La- 
^rangian for the case is obtained by 



£ = q{ilvD v - m)q 
+ G[(qq) 2 + (qi^fq) 2 
~G v (q^q) 2 +U^), 



\q Ci~f 5 T 2 t 2 q\ 



(21) 



where q is the two-flavor quark field, m is the current 
quark mass and t, and Tj are Pauli matrices in color and 
flavor spaces, respectively. In the limit of m = (J, = 0, 
two-color QCD has the Pauli-Gursey symmetry [26, 27], 
so the PNJL Lagrangian is so constructed as to have the 
symmetry. Note that the vector-type four-quark interac- 
tion (<77 M <?) 2 does not work at fi = in the mean-field 
level, since the vector- type condensate is zero at /i = 0. 
The potential U is a function of the Polyakov loop $ and 
the explicit form is shown later in Sec. IV. 

Using the mean-field approximation, one can get the 
thermodynamical potential ft as [28] 



ft = -2N{ 



/||rE[^ ± +^r+>»/ + ) 



with 



1 + 2$e 



-2{iEl 



U = G(a 2 + A 2 ) - G v 



(22) 

(23) 
(24) 



for the chiral condensate a = (qq), the diquark con- 
densate A = \(q 7 ' Ci^T 2 t 2 q)\ and the vector condensate 
(quark number density) n q = (q^q). Here the factors E p 
are defined by 



E* = sgn(E p ± » V )^(E P ± M v) 2 + (2GA) 2 , (25) 
for finite A and 

£± = E p ± fi v (26) 

for A = 0, where sgn(E p ± yu v ) is the sign function, E p = 
M = m — 2Go and /j, v = fi — 2G v n q . In the 
limit of m = fi = 0, the condensate n q is zero, so that 
ft becomes invariant under the rotation in the cr-A plane 
as a consequence of the Pauli-Gursey symmetry. 

In the Polyakov gauge, the Polyakov-loop $ is obtained 
by 



$ = 



i(e* + e-*)=c«#) (27) 



for real number 0, indicating that $ is real. The mean 
fields X = a, A, n q , $ are determined from the stationary 
conditions 



^=0 

ax ' 



(28) 



where ft is regularized by the three-dimensional momen- 
tum cutoff 



d 3 p 
(2tt) 3 



1 

2^2 



dpp 2 



(29) 



because this model is nonrenormalizable. 

At imaginary chemical potential /1 = i9T, ft is invari- 
ant under the extended Z 2 transformation 



+ 7T. 



(30) 
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This can be understood easily by introducing the mod- 
ified Polyakov-loop ^ = e 10 ^ and its conjugate ty* = 
e ~ l6l <j> invariant under the extended Z 2 transformation. 
The condensate A is zero at imaginary /it, since A be- 
comes finite only for fi 2 > m 2 /4 [13, 16]. When A = 0, 
is rewritten into the form of (22) with 



f+ = 1 + 2^e~ 2W ' e -P(Ep-2G v n q ) + e - 



(31) 



f~ = 1 + 2^*e 2W e' p( - Ep+2Gvn "' ) + e~ 2pE ? . (32) 

Obviously, Eqs. (31) and (32) show that fl is invariant 
under the extended Z 2 transformation and at the same 
time that f2 has a periodicity of it in 0, i.e., the RW 
periodicity 



ft(0) = Q(0 + tt). 



(33) 



It should be noted that if A is finite, f2 will not have 
the RW periodicity. Since 0-depcndence of mean fields 
X = <T, Tiq , \& are determined from fi by the stationary 
conditions (28), all the X have the RW periodicity 



X(0) = X(0 + tt). 



(34) 



Furthermore, the RW periodicity ^(9) = ^>(6 + tt) yields 
the relation 



$(0) = -$(0 + tt). 



(35) 



We refer the relation (34) as to Z2-even and the relation 
(35) as to Z2-odd in this paper. 

The PNJL Lagrangian C is invariant under the com- 
bination of the C and the parameter transformation 
—> — 0. This property guarantees that f2(0) = fi(— 0) 
and shows that 9 behaves as a C-odd quantity. Namely, 
C-even quantities such as a, A and $ are 0-even, whereas 
C-odd condensates such as n q are 0-odd. In the PNJL 
Lagrangian, appears only through the pure imaginary 
factor id. This shows that 0-odd quantities become pure 
imaginary, while 0-even quantities become real. The re- 
lation (17) is obtainable from the fact that $(0) is real, 
because 



*(0)* = ($(0)e lS )* = $(-0)e" 



*(-0). (36) 



The PNJL model thus has the same properties as two- 
color QCD for the RW periodicity and the 0-parity. 

The Polyakov-loop $ is 0-even and Z 2 -odd. This indi- 
cates that 



$(0) = -$(tt- 0) 



(37) 



and then $(0)=Oat0 = 7r/2, if $ is a smooth function 
of 0. When T is larger than some temperature T c , ft is 
not a smooth function of at = tt/2. This is the RW 
phase transition. Once the RW phase transition occurs, 
4>(0) is not smooth at = tt/2 and consequently $(0) 
becomes finite there. This means that the spontaneous 
breaking of Z 2 symmetry is induced by the RW phase 
transition. The relation (36), meanwhile, indicates that 



Im^] is 0-odd (C-odd) and then an order parameter of 
the C symmetry. The fact that Im[\E f ] is 0-odd and Z 2 - 
even leads to the same relation 



lm[*(0)] = -Im[*(7T - 



(38) 



as (37). Both the Z 2 and the C symmetry are thus broken 
spontaneously by the RW phase transition, so that the 
two spontaneous breakings take place simultaneously. 



IV. NUMERICAL RESULTS 
A. Parameter setting 

For the N c = Nf = 2 case, we do not have enough 
LQCD data at imaginary chemical potential. We then 
make qualitative analyses here. It is well known from 
the nonlocal version of the PNJL model [5, 6, 29, 30] 
that there is the correlation between the chiral order- 
parameter and the Polyakov-loop in the coupling con- 
stant through the distribution function. This feature 
can be phenomenologically introduced in the local PNJL 
model by using the entanglement vertex [14]. This en- 
tanglement is taken into account in the present analysis. 
We make the following parameter setting. 

1. Since the ratio r = G v /G is of order (iV c )° in 
the leading order of the 1/N C expansion, we take 
r = 0.4 that is determined from the A^ c = 3 case by 
comparing the result of the nonlocal PNJL model 
with LQCD data at finite imaginary chemical po- 
tential. 

2. We introduce the entanglement vertex of the form 
Gi(l — a$ 2 ) for Gi = G and G v , respecting the 
Z 2 symmetry. Here the entanglement parameter a 
is treated as a free parameter, but is assumed to 
be common for both G and G v . We mainly use 
a = 0.4 at which pseudo-critical temperatures of 
the chiral and deconfinement transitions almost co- 
incide when (i = 0. 

3. In the leading order of the 1/N C expansion, is 
scaled by (N c )° and the pion decay constant 

by \/Nc. These scaling properties are assumed to 
determine the parameters G, A and m of the NJL 
sector, where m is simply assumed to be 5.4 MeV. 
In this parameter set, the dynamical quark mass 
M becomes M = 305 MeV at T = fi = 0. The 
resulting parameters are shown in Table I, together 
with the values of m n , and Mo. 

4. Following Ref. [28], we take the Polyakov-loop ef- 
fective potential of the form 



U{<S>) 
T 



= -b 



24e~ a/T $ 2 +m(l - $ 2 



(39) 



with a = 858.1 MeV and b 1 / 3 = 210.5 MeV. This 
potential yields a second-order deconfinement tran- 
sition in the pure gauge limit. 
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The Polyakov-loop potential U used here is determined 
by using the strong coupling expansion of the pure Yang- 
Mills theory. The logarithmic part comes from the Haar 
measure. Since the parameter fitting procedure does 
not refer to microscopic dynamics, it is unclear how the 
Polyakov-loop potential is related to non-perturbativc 
characteristics near T c . This problem should be inves- 
tigated elsewhere by considering other approach based 
on gluon and ghost propagators [31]. 



(a) 

r-0.4 





- (b) 




0.5 






1 ° 






-0.5 







m n [MeV] 


U [MeV] 


M [MeV] 


140 


74.5 


305 


G [GeV"2] 


A [MeV] 


m [MeV] 


7.23 


657 


5.4 



TABLE I: Summary of parameters and physical values. 



B. 9 dependence of order parameters 

Figures l(a)-(c) represent 9 dependence of M, Im^] 
and $, respectively, for three cases of a = 0, 0.2 and 
0.4. Here we consider the high-T case of T = 2.5m,,-. 
Quantities M, lm[9] and $ are order parameters of the 
chiral, the C and the Z2 symmetry, respectively. #-even 
and Z2-odd quantities such as <& are antisymmetric with 
respect to the 9 = tt/2 axis, whereas 9-even and Z2-even 
quantities such as M are symmetric with respect to the 
axis. In addition, 0-odd and Z2-even quantities such as 
Im[^] arc antisymmetric with respect to the 9 = tt/2 axis 
and zero at 9 = and n. As shown in panel (b), Imfvf'] 
has a gap at 9 = tt/2, indicating that the C symmetry 
is spontaneously broken there for the high-T case. The 
zeroth-order discontinuity (gap) in Im^] means that a 
first-order phase transition takes place at 9 = tt/2. This 
is the RW phase transition [7]. At 9 = tt/2, $ also has 
a gap, indicating that the Z2 symmetry is spontaneously 
broken there. The C and Z2 symmetry are always simul- 
taneously broken on the RW phase transition line, since 
Imf*] = sin(6>)$ and hence Im[*] = $ at 9 = tt/2. This 
is a characteristic of two-color QCD. 

The 9-even quantity M has a cusp at 9 = tt/2, when 
the 0-odd quantity Im[4'] has a gap there. In general, the 
zeroth-order discontinuity (gap) in a 9-odd quantity is 
propagated to other 0-odd quantities as the zeroth-order 
discontinuity (gap) and to 0-even quantities as the first- 
order discontinuity (cusp) [32, 33]. The order parameters 
Im[^] and $ are less sensitive to a than M. 



FIG. 1: ^-dependence of order parameters M, Im[\P] and <l>at 
T — 2.5m n . Here M is normalized by the value Mo at T = 
fi — 0. The dotted, dashed and solid lines represent the results 
of q = 0, 0.2 and 0.4, respectively. 

fincmcnt transition, i.e. the Z2-symmctry breaking tran- 
sition, is the second order for a = and 0.2, but becomes 
the first order for a = 0.4. Here the vertical thin-dotted 
line denotes the critical temperature of the first-order 
transition. Since the spontaneous breakings of the Z2 
and C symmetries take place simultaneously, the order of 
the deconfincmcnt transition is the same as that of the 
RW transition at the endpoint. The chiral transition is 
crossover for a = and 0.2, but it becomes the first order 
for a = 0.4. Thus the entanglement intensifies both the 
chiral and deconfinement transitions. For a = and 0.2, 
M has a cusp at T = T| as a result of the propagation 
of the first-order discontinuity in |<I>| [33]. 




T / m, T / m. 



C. Interplay between chiral and deconfinement 
transitions 

Figure 2 shows T-dependence of M/M and |$| at 
9 = tt/2 for three cases of a — 0, 0.2 and 0.4. The decon- 



FIG. 2: T-dependence of (a) M/M and (b) |$| at 6 = tt/2. 
The dotted, dashed and solid lines represent the results of 
a — 0, 0.2 and 0.4, respectively. 

Figure 3 shows T-dcpcndence of M/Mq and $ at 9 = 



6 



(p = 0). Both the deconfinement and chiral transitions 
keep crossover for a = 0, 0.2 and 0.4, although the transi- 
tions become stronger as a increases. The pseudo-critical 
temperature T£ of the deconfinement transition is less 
sensitive to a than the pseudo-critical temperature T° 
of the chiral transition, where the pseudo-critical tem- 
peratures are defined by peak positions of dM/dT and 
rf$/c?T, respectively. 
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FIG. 3: T-dependence of (a) M/M and (b) $ at /i = 0. See 
Fig. 2 for the definition of lines. 



The system at 9 = tt/2 has higher symmetries than 
that at fx = 0. Hence the order parameters have discon- 
tinuities of either the zeroth or the first order at 9 = tt/2. 
The order parameters are considered to be more sensitive 
to a near the discontinuities. This is really true as seen 
by comparing Fig. 2 with Fig. 3. Through the discon- 
tinuities, we can then investigate clearly how strong the 
entanglement between the chiral and deconfinement tran- 
sitions is. The symmetry breakings at the lower bound 
/i 2 = —(Ttt/N c ) 2 thus give deeper understanding. 

Figure 4 shows a-dependence of two (pseudo)critical 
temperatures T£ and T° at 9 = tt/2 and 0. The two 
(pseudo)critical temperatures approach each other as a 
increases, and finally agree with each other at a > 0.2 
for 9 = tt/2 and at a > 0.4 for 9 = 0. As for a < 0.2, the 
speed of the approach is much faster for 9 = tt/2 than 
for 9 = 0. The entanglement thus makes it stronger the 
correlation between the chiral and deconfinement transi- 
tions. Therefore, the difference \T£ — T|| is a good quan- 
tity to determine the value of a. For 9 = tt/2, the decon- 
finement transition is the second order at a < a c sa 0.33, 
but the first order at a > a c . The chiral transition is, 
meanwhile, crossover at a < a c , although it is the first 
order at a > a c . Hence, the point at a = a c is a tri- 
critical point (TCP) for the deconfinement transition and 
a critical endpoint for the chiral transition. 

Figure 5 shows r-dependence of M/Mq and |Im ^| at 
9 = tt/2. The effect of r is similar to that of a, but the 
former effect is smaller than the latter one, when r is 
varied within a realistic range from 0.25 to 0.5. 



FIG. 4: a-dependence of (pseudo)critical temperatures T£ 
and T£ at (a) 9 = tt/2 and (b) [i, = for the case of G v = 0.4G. 
The closed diamond symbols stands for the tri-critical point 
(TCP). The numerical values have ambiguities of about 1 
MeV and the solid lines are drawn by connecting two neigh- 
borhood points. 




FIG. 5: T-dependence of (a) M/M and (b) |Im *| at 9 = 
tt/2. The dotted, dashed and solid lines denote the results for 
r = 0.25, 0.4 and 0.5, respectively. Here a is set to 0.4. 



D. Phase diagram at imaginary and real chemical 
potentials 

First we consider the phase diagram at imaginary [i for 
two cases of a = and 0.4. Here we take r = 0.4. Figure 
6 shows the phase diagram in the 9-T plane. In the left 
panel for the case of a = 0, the upper and lower dotted 
lines denote chiral and deconfinement crossover transi- 
tions, respectively. The two transitions are thus sepa- 
rated from each other, when the correlation between the 
two transitions is weak. In this situation, the C symme- 
try breaking at the RW endpoint is the second order. In 
the right panel for the strong correlation case of a = 0.4, 
the two crossover transitions (dotted lines) agree with 
each other and the C symmetry breaking at the RW end- 
point becomes the first order. In other words, the RW 
endpoint becomes a triple-point where three first-order 
transition lines meet. Thus, the RW endpoint becomes a 
triple-point, when the correlation between the chiral and 
deconfinement transition are strong. 

The Polyakov-loop potential used in this study yields 
the second-order Z 2 symmetry breaking in the heavy 
quark limit. Since the entanglement parameter a makes 
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FIG. 6: The phase diagram in the 0-T plane for a = 
(left panel) and (b) a — 0.4 (right panel). Here the case 
of r = 0.4 is taken. The dotted and solid lines stand for 
crossover and first-order transitions, respectively. In the left 
panel, the upper and lower dotted lines mean the chiral and 
deconfinement crossover lines, respectively. In the right panel, 
the chiral and deconfinement crossover lines almost coincide 
with each other, and the diamond symbol denotes the triple- 
point. 
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the symmetry breaking stronger, the first-order Z2 sym- 
metry breaking in the physical quark mass does not come 
from the Z2 symmetry breaking only. The chiral transi- 
tion becomes stronger, if the coupling constant G is weak- 
ened at finite T; see for example Ref . [34] . In the present 
model, the coupling constant is weakened through the en- 
tanglement vertex, so that the chiral transition becomes 
the first-order. This behavior is seen in the left panel of 
Fig. 2. Therefore, the first-order Z2 symmetry breaking 
is induced by the first-order chiral transition. 

Next we investigate an influence of a on the phase dia- 
gram at real \x. Figure 7 shows T-dependence of M/Mq, 
A/00 and n q for the case of r = 0.4 and /i = m^, 
where <jq is the absolute value of cr at T = /j, = 0. As 
shown in panel (c), A/ctq vanishes around T = 1.5m„-, 
indicating that the superfluid/normal transition occurs 
there. The phase boundary is rather sensitive to a. 

In this study, we take r = G v /G = 0.4 for all the 
calculations. Although this value is obtained from the 
N c = 3 case by comparing the PNJL results with the 
corresponding LQCD data, it is not easy to determine 
the value definitely. We then check r-dependence of order 
parameters, as shown in Fig. 8 where r = 0.25 and 0.5 
are taken as lower and upper limits of a reliable range of 
r. Effects of r on M, $ and A are rather small, although 
r gives an appreciable effect on n q . 

Finally we consider the phase diagram at real /i. Fig- 
ure 9 shows the BCS-BEC crossover (M = /i), the decon- 
finement crossover ($ = 0.5) and the superfluid/normal 
transition lines for two cases of a = and 0.4; see 
Ref. [28] for the definitions of these transitions. When 
T is small, the superfluid/normal transition occurs at 
fx = as expected. When T > and [i > m w /2, 

meanwhile, the superfluid/normal transition takes place 
around T = 1.5m T . The transition line depends on a 
rather strongly compared with the BCS-BEC crossover 



FIG. 7: T-dependence of (a) M/M , (b) (c) A/a and (d) 
n q for the case of r = 0.4 and /j, — m^. The dotted and solid 
lines represent the results of a — and 0.4, respectively. 

and the deconfinement crossover line. Thus the correla- 
tion between the chiral and deconfinement transitions is 
important not only at the lower bound fj 2 = —(Tir/2) 2 
but also at large real /j, such as fi 2 > (m v /2) 2 . 



V. SUMMARY 

Wc have studied properties of two-color QCD at imag- 
inary /j, from the view point of the RW periodicity, the 
charge conjugation and the pseudo-reality. The parame- 
ter n 2 has a lower bound: [i 2 > —(ttT/2) 2 . At the lower 
bound, the system has higher symmetries, that is, the 
Z2 and C symmetries. These symmetries arc simultane- 
ously broken, since their order parameters, $ and Im[4'], 
are identical in virtue of the pseudo-reality The pseudo- 
reality also guarantees that the 1 + 1 system with finite /i 
is identical with the 1 + 1* system with the same amount 
of isospin chemical potential /i, so , so that /i| so has a lower 
bound: iJ? iso > -(ttT/2) 2 . 

The PNJL model has the same properties as two-color 
QCD for the RW periodicity, the charge conjugation and 
the pseudo-reality. The PNJL model is thus a good 
model to investigate two-color QCD at imaginary /i con- 
cretely. We have then investigated the nontrivial correla- 
tion between the deconfinement and chiral transitions at 
imaginary /z for the two-flavor case. The Polyakov loop 
as an order parameter of the deconfinement transition is 
Z2-odd and C-even (0-even) , while the chiral condensate 
as an order parameter of the chiral transition is Z2-even 
and C-even. The order parameters are thus different 8 de- 
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pendence. At 9 = 7r/2, the simultaneous breaking of Z2 
and C symmetries takes place at higher T, and the break- 
ing is induced by the RW phase transition. The order of 
the Z2 and C symmetry breakings at the RW endpoint 
is nontrivial. It cannot be determined by the Z2 and C 
symmetries and the pseudo-reality. The order depends 
on the strength of the entanglement parameter a, i.e., 
the strength of the correlation between the chiral and Z2 
symmetry breakings. The order is the second-order for 
small a, but becomes the first-order for large a. The 
second-order nature is originated in the Polyakov poten- 
tial. Meanwhile, the first-order nature comes from the 
fact that the chiral symmetry breaking becomes the first 
order as a consequence of the strong entanglement. The 
order of the Z2 and C symmetry breakings at the RW end- 
point is thus sensitive to the strength of the correlation 
between the chiral and Z2 symmetry breakings. Finally, 
we have investigated the impact of a on the phase dia- 
gram at real /x. The diagram, particularly the supcrfluid 
transition, is rather sensitive to a. The determination of 
a is thus important for both real and imaginary \x. 



FIG. 8: T-dependence of (a) M/M , (b) $, (c) A/a and (d) 
n q for the case of a — 0.4 and [i — m^. The dotted, solid and 
dashed lines represent the results of r = 0.25, 0.4 and 0.5, 
respectively. 




1 2 



r = 0.4 

a = 0.4 
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At the present stage, we do not know how large a is, 
but it is possible to determine the value of a from LQCD 
simulations at 9 = 7r/2, particularly by seeing the correla- 
tion between the chiral and Z2 symmetry breakings. It is 
interesting as a future work. In two-color QCD, the cor- 
relation between the chiral and deconfinement transitions 
is thus important for both real and imaginary /i. This 
is true also for three-color QCD. This strongly suggests 
that understanding of three-color QCD at imaginary [i is 
important to determine the phase diagram at real [i. 



FIG. 9: The phase diagram in the fi-T plane for a = (left 
panel) and 0.4 (right panel). Here the case of r = 0.4 is taken. 
The dotted, dashed and solid lines stand for the BCS-BEC 
crossover (M = /j,), the deconfinement crossover ($ = 0.5) 
and the sup erfluid /normal transition lines, respectively. 
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